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ABSTRACT 

We constrain the uncertainty in waiting times for detecting the first double-neutron- 
star (DNS) mergers by gravitational wave observatories. By accounting for the Poisson 
fluctuations in the rate density of DNS mergers and galaxy space density inhomogene- 
ity in the local Universe, we define a detection 'zone' as a region in a parameter space 
constrained by the double neutron star merger rate and two LIGO operations param- 
eters: an observation horizon distance and science run duration. Assuming a mean 
rate of about 80 DNS mergers per Milky Way galaxy Myr -1 , we find a 1/20 chance 
of observing a merger by Enhanced LIGO in only 1 yr of observation. The minimum 
waiting time and temporal zone width for an Advanced LIGO sensitivity are much 
shorter and imply that there is a 95% probability of detecting a DNS merger in less 
than 60 days and a 1/20 chance of a first detection in about 1 day. At the 5% probabil- 
ity threshold for a first detection, we find that the effect of galaxy clusters on detection 
is smoothed out and may only influence detection rates after 5-10 years observation 
time. 

Key words: gravitational waves - binaries: neutron stars - gamma-rays: bursts 



1 INTRODUCTION 

Double neutron star (DNS) binary mergers are potentially 
strong sources of detectable gravitational wave (GW) emis- 
sion. Nine DNS binary systems have been discovered, eight 
of them are in the Galactic disk and one is in a globular 
cluster (Stairs 2004, Lorimer et al. 200 6). Energy loss from 
GW emission (|Tavlor fc Weisberel 19891 ) causes an orbital in- 
spiral until the binary system merges, resulting in a burst of 
GWs usually described as a 'chirp' signal. 

The US LIGO (Laser Interferometer Gravitational-wave 
Observatory) is searching for, among other potential sources, 
DNS mergers in the local Universe and has extended its sen- 
sitivity horizon to these events to some tens of Mpc (Abbott 
et al. 2007). LIGO sensitivity to a DNS merger depends on 
its distance and relative orientation to the event. To aid 
searches a catalogue has recently been compiled of distances, 
sky positions and blue luminosities of galaxies in the local 
Universe: t he Compact Binary Co alescence Galaxy (CBCG) 
catalogue (|Kopparapu et al]|2008l ). 

Theoretical calculations and simulations predict a wide 
range of delay times between DNS formation and merger 
(see Belczynski et al. 2006). Their results support the ar- 
gument that a significant fraction of DNS mergers follow 
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massive star formation. Blue luminosity is a known indica- 
tor of massive star formation and may also provide a means 
of tracking the DNS merger rate, for binaries with merger 
times that are short compared to the Hubble time. The rel- 
atively frequent occurrence of such systems is supported by 
the discovery of the double pulsar system J0737-3039A,B. 
With an orbital period of only 2.45 hr, it will coalesce in 
only 87 Myr. In the local Universe, out to some hundreds 
of Mp c, we assume, like most other studies, e.g. IPhinnevI 
(1991) and others, that the DNS merger rate scales with 
blue luminosity. 

Interestingly, some short gamma ray bursts (GRBs), ob- 
served by the Swift satellite, may also be the high-energy 
EM emissions from DNS mergers. This is evidenced by the 
optical localization of several short hard bursts with their 
host galaxies. The identification of short GRBs with differ- 
ent types of galaxies and distances from galaxy centers is 
related, in a non-trivial way, to the DNS velocity and the 
delay times between formation and merger (see e.g. Belczyn- 
ski et al. 2006) . In contrast to the binary pulsar systems, the 
offset position for GRB050509b at 40 ±13 kp c from its ellip- 
tical galaxy G17 implies large kick velocities (|Grindlav et al.l 
2006). This discrepancy between the observed radio pulsar 
binaries and short GRB locations could indicate that DNS 
mergers may occur in both old stellar populations, such as 
globular clusters, and relatively younger environments. Even 
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though the exact distribution of DNS mergers with galaxy 
types is not well constrained, it is reasonable to assume that 
there is a strong correlation between massive star forma- 
tion rates, DNS merger rate and the blue luminosity of the 
progenitors. 

With the assumption that galactic blue luminosity is a 
direct measure of massive star formation in the local Uni- 
verse, the CBCG catalogue can be used to infer the DNS 
merger rate at extra-galactic distances by scaling from the 
blue luminosity of our Galaxy. This requires knowledge of 
the Galactic DNS merger rate. In order to estimate the 
DNS coalescence rate, Kalogera et al. (2001) used a semi- 
empirical approach, based on the observed properties of 
known DNS systems and pulsar survey selection effects, to 
obtain scale factors that correct for the unobserved fraction 
of existing systems. With this model and more recent pul- 
sar su rvey results that include J0737-3039, iKalogera et al.l 
(2004) presented bounds for the merger rate for Galac- 
tic disk DNS sy s tems, yielding Hvns = SStf^ 1 Myr -1 . 
iRegimbau et al.l (|2005T ). using numerical simulations, find 
that mergers are possible between 2 x 10 5 yr and the age 
of the Universe. Using both evolutionary and statistical 
models, IRegimbau et al l (j2005T l find a Galactic merger rate 
of 17 Myr -1 , s imilar to the lower bound calculated by 
IKalogera et"atl l|2004h . It is clear from these studies that 
the DNS merger rate is highly uncertain by several orders 
of magnitude. To account for this uncertainty, we employ a 
lower range of values for the DNS merger rate of 1 — 100 
Myr -1 per Milky Way Galaxy. 



2 THE MODEL 

2.1 DNS merger detection rates 

In a previous study, ICoward et al.l l|2005t ) employed a sim- 
ple Euclidean model for the rate density of DNS mergers to 
study the fluctuations in the detection rate of DNS by LIGO 
and Advanced LIGO. In this study, the aim is to show how 
the waiting time for a first detection of a DNS merger by 
next generation LIGO depends on the distribution of galax- 
ies, the sensitivity of the detector and the observation time. 
To model the effect of the uncertainty in waiting times of 
DNS mergers, we start with a semi-empirical model to ac- 
count for fluctuations in the cumulative merger rate arising 
from inhomoge neities in the local Univer se. 

Following iKopparapu et all {2008), we employ the 
CBCG catalogue, using the cumulative sum of blue lumi- 
nosities from the catalogue, Cl{D), where D is distance. 
We re-scale this galactic blue luminosity as a function of dis- 
tance to detector horizon distance Dh,LJso as to account for 
the response of the LIGO detector to GWs from a randomly 
oriented DNS merger averaged over all sky positions. The cu- 
mulative number of DNS mergers, A r DNs(-Dh), in observation 
time T detectable within horizon distance < 100 Mpc is 
approximated using Cl(D^). For Dh > 100 Mpc, the galaxy 

1 This distance is denned as the maximum distance at which 
LIGO can detect an optimally orientated DNS merger with signal- 
to-noise ratio of eight. For non-optimally orientated mergers, the 
sensitivity distance, D s , accounts for the antenna pattern of the 
detector, and is denned as D s = D^/\fE 
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Figure 1. Cumulative DNS mergers detected as a function of 
LIGO horizon distance assuming the most probable rate density 
of 83 mergers per Milky Way galaxy yr -1 from Kalogera et al. 
(2004) and the scaling of blue luminosity to DNS number density 
from Kopparapu et al. (2008) who used the Compact Binary Co- 
alesence Galaxy catalogue. We extrapolate the cumulative num- 
ber beyond 100 Mpc using a simple Euclidean cubic power law 
(dashed curve). 

distribution is considered uniform, so we extrapolate using 
a cub ic power law, similar to eq. (9) in IKopparapu et al.l 
|200Sl ). Thus, for D h < 100 Mpc, 

™— -(z#^)(^) (£).<>> 

and, for D h > 100 Mpc, 

where Lio = 10 10 Lb, q is in units of solar blue luminosity. 

Figure [1] plots equations 1 and 2 assuming an obser- 
vation time of 1 yr and a Milky Way DNS merger rate of 
83 Myr -1 , correspon ding to the most probable rate from 
IKalogera et alj (|2004T ). For horizon distances comparable to 
that of Enhanced LIGO, ~ 60 Mpc, the rate is about 1 
per 20 years; for a horizon distance comparable to that of 
Advanced LIGO, ~ 450 Mpc, the rate is about 20 per yr. 



2.2 Waiting times for detection of DNS mergers 

In addition to the uncertainty in the mean DNS merger rate, 
detection also depends on the intrinsic uncertaint y in the 
waiting times between events. iHowell et all l|2007r i demon- 
strated, using GRB peak flux data from the Swift satellite, 
that the waiting times of cosmic transients follow Poisson 
statistics. The study derived a power-law relation that is 
independent of the luminosity distribution of the sources to 
predict waiting times for bright GRBs observed by Swift. At 
the core of these models is the inclusion of the Poisson-based 
uncertainty in observation waiting times for astronomical 
transients. 
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Table 1. The waiting times for a first detection, in days, for 
the 5% and 95% probability thresholds assuming horizon and 
sensitivity distances, and D s , corresponding to Advanced and 
Enhanced LIGO, along with the temporal zone widths. 





T 5% 


T 95% 


Tzonc 


Adv. LIGO D h 


1 


60 


60 


Adv. LIGO D s 


10 


600 


590 


Enh. LIGO D h 


365 


21900 


21800 


Enh. LIGO D s 


2000 


< 10 4 


< 10 4 



As in that work, we assume that DNS mergers are inde- 
pendent events, i.e. the occurrence of a DNS merger in one 
region of space does not influence the probability of another 
such event occurring. Therefore, the temporal separation be- 
tween events detected by LIGO will follow an exponential 
distribution, defined by a mean event rate above 7?.DNs(£ ) h) 
for events out to some horizon distance Dh in Euclidean 
space. The probability for at least one event to occur in a 
volume bounded by Dh during an observation time T at 
constant probability e is given by: 

V(n > 1;7£dns(Ai),T) = 1 - e KDBs(Di)T = e. (3) 

For this equation to remain satisfied with increasing obser- 
vation time: 

ATdnsCAO = R DN s(D h )T = |ln(l - e)| , (4) 

with mean number of events Nr>Ns{Dh) ■ 

The probably event horizon (PEH) is defined by the 
minimum detectability horizon distance, Df EH (T), for at 
least one event to occur over some observation time T, with 
probability above some selected threshold e. L»f BH (T) is de- 
fined by by fixing e and solving the above relation, thus 
defining the PEH l|Coward et al.ll2005l ). 

In this study, we use two thresholds e = (0.5,0.95), 
corresponding to 5% and 95% probability respectively of 
observing at least one event within D^. In other words we 
would be surprised to observe an event before T(D^% n ), 
henceforth defined as the 'waiting time', and extremely sur- 
prised if no events were observed by T(Dq^). For obvious 
reasons we denote the interval between these two times the 
'detection zone': 

WAO = T(D^) - T(D^ H ) . (5) 



3 RESULTS AND DISCUSSION 

Figure 2 plots the detection zone defined by equation ©. 
The Dg^ H curve shows that the inhomogeneity of the lo- 
cal Universe affects the PEH at horizon distances less than 
50 Mpc. Waiting times corresponding to < 50 Mpc 
are reduced compared with a uniform galaxy distribution 
(dot-dashed line) but the deviation manifests after rela- 
tively long observation times i.e. 5—10 yr. For a horizon 
distance corresponding to a sensitivity of Enhanced LIGO, 
~ 60 Mpc, the waiting time is about 1 yr and the detec- 
tion zone spans about 60 yr. Assuming a horizon distance 
of Advanced LIGO, the waiting time is about 1 day and the 
detection zone spans 60 days. This implies that there is a 
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Figure 2. The plot shows the evolution of the DNS merger detec- 
tion zone (shaded region) bounded by the 5% and 95% probability 
contours for detection by LIGO. The intersection of a horizon dis- 
tance (vertical axis) with the 5% probability contour defines the 
waiting time before entering the detection zone. Assuming the 
most probable rate density of about 80 DNS mergers per Milky 
Way galaxy yr - 1 from Kalogera et al. (2004), the minimum length 
of the science runs to enter the zone is 365 and 2 days for En- 
hanced and Advanced LIGO sensitivities respectively. The detec- 
tion zone width for Advanced LIGO is constrained to 2—60 days, 
while for Enhanced LIGO it is 1—60 yrs. The deviations of the 
probability contours from a simple power law, dot-dash curve, at 
horizon distances less than 50 Mpc is a result of fluctuations in 
the space density of galaxies and clustering in the local Universe 
(see Figure 3). 




5 10 15 20 

observation time (yr) 

Figure 3. Plot of the fractional difference between the 5% detec- 
tion probability contours (PEH 5 %) assuming a uniform distribu- 
tion of galaxies compared to a model that includes fluctuations in 
the space density of galaxies in the local Universe. Clustering in 
the local Universe results in a maximum decrease of about 1.4 in 
the PEH 5 % distance at an observation time of about 7-10 yr. For 
practical waiting times for a first detection i.e. 1-5 years, the effect 
from inhomogeneities of the galaxy distribution is not significant. 
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Figure 4. Plot of the detection zone, assuming an Advanced 
LIGO horizon distance of 450 Mpc, with Galactic DNS merger 
rates (vertical axis) spanning 1 - 100 Myr -1 per Milky Way 
galaxy. The intersection of a particular rate with the 5% and 95% 
probability contours gives the waiting times. The corresponding 
waiting times for the upper and lower rates at the 5% level are 
0.5 and 90 days respectively. 



90% probability of observing at least one DNS merger in 
1 — 60 days of a science run at Advanced LIGO sensitivities. 

Figure 3 plots the fractional difference between the 
PEH 5 % contours for a model that assumes a uniform galaxy 
distribution and one that includes inhomogeneities in galaxy 
density in the local Universe. There is a maximum factor of 
about 1.4 decrease in the horizon distance from clustering in 
the local Universe, corresponding to an observation time of 
7-10 years. Considering that Enhanced and Advanced LIGO 
should have horizon distances out to the nearly homogenous 
galaxy space density regime, it is clear that this effect will 
not be important for science runs of practical duration. 

iKopparapu et all l|200ct ) show that the relatively small 
horizon distances of the two initial LIGO detectors can vary 
because of the relative orientation of the detector to the sky 
positions of galaxies in the local Universe. This is because 
the detector antenna pattern must coincide with the sky po- 
sition of the DNS host galaxies for a detection. The galaxies 
become sparser and the sky separations become larger as 
the horizon distance becomes smaller. In this analysis, the 
effect of different detector locations has not been included, as 
this study focusses more on detector sensitivities comparable 
to those of Enhanced and Advanced LIGO. The associated 
horizon distances of these detectors extend from the border 
of inhomogenous to mostly uniform galaxy space densities. 

Figure 4 plots the detection zone at a fixed horizon dis- 
tance corresponding to that of Advanced LIGO, but with a 
varying mean DNS merger rate of 1 — 100 Myr -1 per Milky 
Way Galaxy to account for the uncertainty in this rate. The 
corresponding minimum waiting times at the 5% probability 
level for the two rate extremes are 0.5 and 90 days. 

In summary, inclusion of the Poisson uncertainty for 
the minimum waiting times of DNS mergers detectable by 
GW observatories provides a more realistic means for defin- 
ing a detection zone in terms of the duration of a science 



run. Even if the rate density of DNS mergers is well con- 
strained, the study demonstrates how the waiting times can 
vary widely over an observation period. To highlight this, the 
model calculates a 1/20 chance of observing a DNS merger 
by Enhanced LIGO in only 1 yr of observation assuming an 
event rate of about 80 Myr -1 in the Milky Way Galaxy. The 
horizon distance of Enhanced LIGO is also in the distance 
regime where inhomogeneities in galaxy space density will 
only marginally influence the probability of a detection. 

The waiting times and zone width at an Advanced 
LIGO sensitivity are much shorter in duration and imply 
a first detection at the 95% probability level will occur in 
less than 60 days and there is a 1/20 chance of a first detec- 
tion in about 1 day. If the merger rate is of order 1 Myr -1 
in the Milky Way Galaxy, the waiting times at a 5% prob- 
ability level assuming an Advanced LIGO horizon distance 
are extended to about 90 days. 



ACKNOWLEDGMENTS 

D. M. Coward is supported by an Australian Research Coun- 
cil grant (LP0667494). This research was made possible by 
a visit to the Observatoire de la Cote d'Azur, funded by 
the Australian Academy of Science: Scientific visits to Eu- 
rope program. The host institute provided a fertile environ- 
ment for discussions with Dr Tania Regimbau that stimu- 
lated this work. The author also thanks Dr Ron Burman, 
Prof. David Blair for fruitful discussions and Prof. Richard 
O'Shaughnessy for providing a review on behalf of the LIGO 
Scientific Community. Finally, the author thanks the anony- 
mous referee for providing very informative comments and 
suggestions for improving this paper. 



REFERENCES 

Abbott B., et al. 2007, preprint (gr-qc/0704.3368) 
Belczynski K., Perna R., Bulik T., Kalogera V., Ivanova 

N., Lamb D.Q., ApJ, 2006, 648, 1110 
Coward D.M., Burman R.R. MNRAS, 2005, 361, 362 
Coward D.M., Lilley M., Howell E.J., Burman R.R, Blair 

D.G., 2005, MNRAS, 364, 807 
Grindlay J., Zwart S.P., McMillan S., Nat. Physics, 2006, 

2, 116, doi: 10.1038/nphys214 
Hulse R.A., Taylor J.H. 1975, ApJ, 195, L51 
Howell E., Coward D., Burman R., Blair D., 2007, ApJ 666, 

L65 

Kalogera V., Narayan R., Spergel D.N., Taylor J.H. 2001, 

ApJ, 556, 340 
Kalogera V., et al. 2004, 614, 137 
Kopparapu R.K. et al., 2008, ApJ, 675, 1459 
Lorimer D.R., Stairs I.H., Freire P.C., Cordes J.M., Camilo 

F., Faulkner A.J., Lyne A. J., Nice D.J. et al., 2006, ApJ, 

640, 428 

Phinney E. S. 1991, ApJ, 380, L17 

Regimbau T., de Freitas Pacheco J. A., Spallicci A., Vincent 

S. 2005, Class. Quantum Grav., 22, S935 
Taylor J.H., Weisberg, J.M. 1989, ApJ, 345, 434 
Stairs I.H., 2004, Science, 304, 547 
Wolszczan A. 1991, Nature, 350, 688 



